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ABSTRACT 

This paper presents a new technique for reconstructing the spatial distributions of hydrogen, tem- 
perature and metal abundance of a galaxy cluster. These quantities are worked out from the X-ray 
spectrum, modeled starting from few analytical functions describing their spatial distributions. These 
functions depend upon some parameters, determined by fitting the model to the observed spectrum. 
We have implemented this technique as a new model in the XSPEC software analysis package. We 
describe the details of the method, and apply it to work out the structure of the cluster A1795. 
We combine the observation of three satellites, exploiting the high spatial resolution of Chandra for 
the cluster core, the wide collecting area of XMM-Newton for the intermediate regions and the large 
field of view of Beppo-SAX for the outer regions. We also make a threefold test on the validity and 
precision of our method by i) comparing its results with those from a geometrical deprojection, ii) 
examining the spectral residuals at different radii of the cluster and in) reprojecting the unfolded 
profiles and comparing them directly to the measured quantities. Our analytical method yields the 
parameters defining the spatial functions directly from the spectra. Their explicit knowledge allows a 
straightforward derivation of other indirect physical quantities like the gravitating mass, as well as a 
fast and easy estimate of the profiles uncertainties. 

Subject headings: galaxies: clusters: general cooling flows — intergalactic medium — galaxies: clus- 
ters: individual (A1795) — methods: data analysis — X-rays: galaxies: clusters 



1. INTRODUCTION 

A deep insight into the early life of our Universe is 
provided by the largest and youngest virialized struc- 
tures: clusters of galaxies. A valuable tool to fathom 
the physics and the structure of these objects is their 
X-ray emission. This stems from the hot (1 — 10 keV) 
and optically thin plasma (also known as intra-cluster 
medium, or ICM) which makes up some 10 — 20% of 
their total mass. This rarefied plasma is so hot because 
it has to develop a pressure strong enough to maintain 
itself in quasi-hydrostatic equilibrium within the gravita- 
tional well of the dark matter. The ICM is transparent 
and allows to see deeply, down to the very core of the 
cluster. This transparency has also a drawback: photons 
born at different depths mix up along the same lines of 
sight, and a direct use of the spectral information to in- 
vestigate in three dimensions the physical conditions of 
a cluster is impossible. The first step is necessarily a 
deprojection of the X-ray spectrum, which provides the 
deconvolved three-dimensional profiles of density, tem- 
perature and metal abundance. From these fundamental 
direct quantities others may follow, like the distribution 
of heavy elements or the cluster gravitating mass. 

This paper describes a novel deprojection technique 
we have implemented as a new model in the XSPEC 
software package for the analysis of X-ray spectra. The 
fundamentals of the theory are recalled in Section [3 In 
Section |3| we apply them to build up our deprojection 
model. In Section 0] we make a specific example by de- 



projecting the nearby cluster A1795. A comparison is 
also made between our results and those from another 
deprojection technique. In Section[5lthe methods to test 
the solidity of our deprojection are reviewed. Finally, in 
Section we show how the result of our analytical de- 
projection may help to derive other physical variables of 
cosmological interest. 

2. THE SPECTRUM OF AN OPTICALLY THIN SPHERE OF 
PLASMA 

In this Section we briefiy review the fundamental for- 
mulae for the flux from the diffuse medium of a cluster 
of galaxies. This is approximated as a spherically sym- 
metric cloud of optically thin plasma located at redshift 
z. The observed photon flux from a source subtending 
the solid angle fl is 



(1) 



Te^ Ie/E, 



where Ie is the specific brightness in the observer's rest 
frame, expressed in terms of the spectral energy E. The 
specific brightness is obtained by integrating the spectral 
emissivity j^/ along an optical path across the source. If 
this passes at a distance b from the emission center, the 
hypothesis of spherical symmetry yields 
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dr- 



• JE' {r) 



\fr^ 



(2) 



1 lASF/CNR Sez, 
Milano - Italy; fabio' 



"G. Occhialini" Via E. Bassini 15/A 20133 
mi.iasf.cnr.it. 

^ Dipartimento di Scienze, Universita dell'Insubria, Via Valleg- 
gio 11, 22100, Como - Italy. 

^ Osservatorio Astronomico di Brera, via Bianchi 46, 22055 Mer- 
ate (LC), Italy. 



where all primed quantities refer to the cluster rest frame, 
so E = _E'/(1 + z) is the redshifted spectral energy. 
The quantity Ie/E^ is invariant along the path of a 
light ray (see e.g. Rybicki & Lightman 1979), thus the 
emitted and observed specific brightnesses are related 
by I'p^E'^ = Ie/E^. Taking this into account. Equa- 
tion Q can be substituted into J^l, where the integration 
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over the solid angle can be usefully rewritten in terms of 
the projected size of the source. If the emitting region 
is located at the angular distance Da and is seen as an 
annulus bounded by the projected rims 6niin and 6niax, 
Equation becomes 

where = js/E is the photon emissivity. The spectral 
emissivity is jE = nn n^-^-E (T, Z) , where and hh are 
the electron and hydrogen densities: the temperature T 
and the metal abundance Z enter jE only through the 
spectral cooling function A^;. The occurrence of a dou- 
ble integral in the evaluation of a spectrum is a compu- 
tationally tough problem. The possibility of writing a 
reasonably efficient XSPEC code to model the spectrum 
of a cluster lies in the possibility of reducing this dou- 
ble integral to a single integral. Under the hypothesis 
of spherical symmetry this can be done, and the reader 
is referred to Appendix A for the technical details. The 
resulting formula 
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is the basis of our deprojection method. 

3. THE PROBLEM OF DEPROJECTION 

3.1. The geometrical deprojection 

The problem of recovering the 3-D shape of T, and 
Z is usually tackled with the methods of geometrical de- 
projection. Our term "geometrical" encompasses a vari- 
ety of non-parametric techniques to extract 3-D quanti- 
ties from their corresponding 2-D projected counterparts. 
This Subsection provides a summary of the main geomet- 
rical deprojection techniques, as well as some references 
for further readings on the subject. 

The method of deprojection in the study of galaxy clus- 
ters was introduced by Fabian et al. (1981), and Kriss, 
Cioffi & Canizares (1983) in their pioneering works on 
the analysis of Einstein data. The spectral resolution 
of the instruments HRI and IPC on board Einstein was 
too poor for a direct deprojection of the spectra, and 
these authors unfolded the cluster starting from its sur- 
face brightness. In the hypothesis of spherical symmetry, 
the image of the cluster is divided into N rings of radii 
ri < r2 < ■ ■ ■ < rN centered on the emission peak. The 
3-D emission is supposed to come from N homogeneous 
spherical shells bounded by the same radii. The observed 
surface brightness of the j-th projected ring (of area Aj) 
is 



of the instrument. Since the emission rate depends both 
on temperature and density, in the absence of further 
spectral information this degeneracy hinders any deter- 
mination of these quantities. In oder to determine sepa- 
rately the 3-D density and temperature profiles, Fabian 
(1981) was forced to assume an average metallicity of 
the cluster, and to suppose the gas to be in hydrostatic 
equilibrium within an assumed gravitational potential. 

When spectral information became available, the hy- 
pothesis of hydrostatic equilibrium could be avoided, 
thus allowing a substantial extension of this deprojection 
technique. Nulsen & Bohringer (1995) adopted a variant 
of Fabian's method to analyze ROSAT spectra: it was 
applied to the count rate of each energy channel, instead 
of the integrated spectrum. This procedure yielded the 
pulse height spectrum of each spherical shell, which was 
converted to photons and fitted to an emission model 
to obtain the 3-D values of temperature and hydrogen 
density. 

David et al. (2001) ~ see also Allen, Ettori & Fabian 
(2001) and Buote (2000) - modeled the fiux from any 
annulus as a weighted sum of the emissivities of the 
spherical shells. Each shell was supposed to be character- 
ized by its own uniform temperature, density and metal 
abundance, considered as the parameters of the fit, de- 
termined by fitting simultaneously the annular spectra. 
™Eltori et al. (2002) devised a matrix method to depro- 
q^fspectra. The spectrum of each ring is first analyzed 
to determine its emission weighted projected tempera- 
ture, metal abundance, emission integral and luminosity 
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where Vij is the volume of the fraction of the i-th spher- 
ical shell subtended by the j-th projected ring. This ma- 
trix equation can be solved for the spatial count emission 
rates Ci of each shell. In turn, these quantities are con- 
verted to photon emission rates with the effective area 



where ~ 1.21nH is the electron density, Vij has the 
same meaning as before, and is the emissivity (in erg 
cm~^ s~^) of the z-th shell. Since the volumes are known, 
the previous matrix equations can be solved to yield the 
3-D quantities appearing on their right-hand sides. 

3.2. The analytical deprojection 

Whenever a spectral information is available, the de- 
projection techniques described in the last Section do not 
require the introduction of any "theoretical bias" in the 
3-D profiles. Nevertheless, the profiles recovered are of- 
ten rather wiggly. Thus, for many applications they need 
to be regularized, causing the loss of some information on 
the fine-grained patterns of the X-ray emitting gas. An- 
other common requirement of all geometrical techniques 
is the lack of gaps between nearby rings, since they would 
bar from knowing the exact number of photon counts to 
subtract and proceed inwards. 

In this paper we introduce an alternative approach to 
deprojection in which this restriction is not necessary. 
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and where no a posteriori regularization is required. The 
only condition our method does not relax is the hypoth- 
esis of spherical symmetry, under which all the relevant 
deprojection formulae discussed in this paper have been 
deduced. Our method looses some fine-grained infor- 
mation (exactly as the geometrical techniques do), the 
difference being that here the loss occurs before the de- 
projection. In our analytical deprojection we model the 
X-ray emission from the cluster, i.e. we do introduce a 
"theoretical bias" by giving to the 3-D distributions of T, 
uh and Z known functional forms, dependent on some 
parameters ai, 02, • • • oat. With Equation we calcu- 
late the model spectrum TE[ai, 0,2, - ■ ■ un] as a function 
of these parameters, fit it to the observed spectra and 
determine the best- fit values of ai, 02, • ■ ■ ajv, fixing the 
deprojected 3-D distributions of T, njj and Z. These 
profiles do not need any regularization, simply because 
-strictly speaking-ours is not a "deprojection", but a 
procedure to determine the parameters of smooth 3-D 
functions which better match the 2-D observed profiles. 
One might question the strength of the bias introduced 
by forcing the 3-D distributions to be described by a func- 
tional form with a (necessarily) small number of param- 
eters. This important issue shall be discussed in Section 

El . 

We have implemented our deprojection technique as 
a new model for the XSPEC analysis software package 
and have christened it smaug after the dragon of J.R.R. 
Tolkien's tale The Hobbit, but the name also stands for 
Spectral Modeling And Unfolding of Galaxy clusters. 
Given a first guess of the fit parameters ai, a2, ■ • • , 
XSPEC calls smaug and calculates for it the spectral 
cooling function Ae{T, Z). Then smaug sums up the lo- 
cal contributions to the spectrum by evaluating the right- 
hand side of Equation (0J, and passes back the result to 
XSPEC. Finally, XSPEC compares the model spectra to 
the observed spectra and works out an improved esti- 
mate of ai, a2, • • • Oat. The procedure is iterated until 
the values of the fit parameters are stable. 

The evaluation of the smaug model requires the cal- 
culation of the integral 10}, deduced from the double 
integral Q under the hypothesis of spherical symme- 
try. This reduction is necessary, since double integrals 
are not easy to evaluate numerically. Besides, the inte- 
grand function contains the spectrum e£;(r), i.e. a func- 
tion not only of r, but also of the energy E. Without 
the reduction of the double integral, the evaluation of 
a model spectrum would be exceedingly time-consuming 
from the computational viewpoint. Even the numerical 
calculation of the single integral requires some care in 
order to keep as light as possible the computational effort 
without loosing accuracy. We have confined all technical 
considerations about this important topic to Appendix 
B, and the interested reader is referred there. 

The model is fed with a collection of spectra extracted 
in annuli centered about the X-ray emission peak. Hav- 
ing specified the cosmology, the user is prompted for a 
series of parameters. The first are the redshift z of the 
source, the rims bmin and bmax (in arcmin) of each extrac- 
tion ring, as well as the cutoff radius (in Mpc) and the 
number of mesh points for the calculation of the integral 
(0}. Two further fixed parameters specify whether the 
spectrum is to be calculated or interpolated form existing 
tables, and which type of plasma emission code has to be 



used: the choice is among Raymond-Smith, Mekal, Meka, 
or APEC. The variable fit parameters of the model define 
the distributions of hydrogen density, electron tempera- 
ture and metal abundance within the cluster. 

The hydrogen density profile is modeled with a double- 
beta function, fixed by six parameters: 



nnir) ^no[fnc{r) + (I-/)%(r)] 



where 



1 



-/3c, B 



(6) 



(7) 



The first beta-component refers to the ICM, and the sec- 
ond to the dominant galaxy (if any) lying at the cen- 
ter of the cluster. The variable uq is the central hy- 
drogen density, and / is the fraction of uq relative to 
the ICM. The rationale for the choice © is that double 
beta-models have been widely and successfully applied to 
fit the density distributions in several clusters (see e.g. 
Mohr, Mathiesen & Evrard 1999). 

The next group of eight parameters defines the tem- 
perature profile 



T(r) = fc(r) [To+Til{r)] 
k{r) = [l + (r/rtaii)2]-'" 

J, _ 2 arctan[(r/nao)''] 
; — TT l+(rcool/i-)« ■ 



(8) 



Despite its appalling appearance, this formula has a sim- 
ple astrophysical rationale. There are three regions with 
characteristic radii riso, t'cooI and rtaii, the hierarchy of 
which is expected to be rjso < TcooI < '"tail- At small 
radii, i.e. for r < riso, the plasma is nearly isothermal, 
being T ^ Tq. This core flattening of T has been re- 
cently observed with Chandra in some clusters (A2199: 
Johnstone et al. 2002; A1835: Schmidt, Allen & Fabian 
2001; A1795: Ettori et al. 2002). At intermediate radii 
Hso S ^ S ''cool the temperature T To + Til{r) in- 
creases with radius, as observed in practically all relaxed 
clusters harboring a cD galaxy at their center (see e.g. 
De Grandi & Molendi 2002). At radii larger than TcooI 
the temperature flattens, and finally for r > rtaii it starts 
to decrease, behaving essentially as a beta function, in 
the same fashion as the hydrogen density. In the out- 
skirts of the cluster T cx r~^" and uh oc r~^^^, thus the 
equation of state of the ICM is expected to be a poly- 
trope T cx n^J^" . The decrease of T at large radii was 
first observed by Markevitch et al. (1998) with ASCA in 
a sample of 30 clusters, and more recently by De Grandi 
& Molendi (2002) with Beppo-SAX. 

It is worth to remark that an apparently complicated 
profile like © is actually highly modular: by suitably 
fixing some parameters, the user may impose the tem- 
perature profile to be only increasing, only decreasing, 
or even isothermal. A function like © is able to encom- 
pass a wide range of qualitative behaviors of the spatial 
temperature profile. 

The metal abundance profile is described by the three 
parameters function 



Z[r) = Zo 



1 



-C 



(9) 



This choice is motivated by the observation of a ra- 
dial decrease in the metal concentration in practically 
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all clusters (De Grandi & Molendi 2001), and by more 
recent Chandra and XMM-Newton observations showing 
a flattening of the metallicity in the core of M87/ Virgo 
(Molendi & Gastaldello 2001) and other clusters (e.g. 
A496: Tamura et al. 2001). For the time being we have 
not included in the functional form of Z{r) the decrease 
in metallicity observed in the core of Perseus (Schmidt, 
Fabian & Sanders 2002) and of Centaurus (Sanders & 
Fabian 2002). If called for by future observation, it will 
be easy to modify the profile (jS)) to properly describe this 
behavior. 

It is known that XSPEC rescales any additive model 
by multiplying it for a normalization constant A/" (Ar- 
naud & Dorman 2000), thus the model parameter rtg in 
Equation lO is degenerate with respect to M . Before 
starting a fit, no must be fixed to unity. The normaliza- 
tion M imposed to smaug by XSPEC has the meaning of 
the square of the central hydrogen density in cgs units, 
so AA^/^ — nnir ^ 0) cm^^. 

As a final remark, we would like to stress that our 
model is extremely fiexible: the code has been written 
with a modular structure, allowing an easy change of 
the functional forms of uh, T ov Z, whenever new and 
better X-ray observations will require it. For instance, to 
modify the temperature profile , the user has to 1) re- 
place the code fragment returning the temperature with 
a more suitable one, 2) edit the parameters input section 
in the main program, and finally 3) edit the Imodel . dat 
ASCII file telling XSPEC what are the model parame- 
ters (Arnaud & Dorman 2000). We hope that the good 
flexibility might render smaug (and its possible future 
implementations) a useful tool for the whole community. 

4. A SPECIFIC EXAMPLE: THE DEPROJECTION OF ABELL 

1795 

4.1. Data preparation 

In this Section we describe how our deprojection code 
works in practice. We deproject the spectrum of A1795, 
a nearby (z = 0.0631) and relatively relaxed cool-core 
cluster (Ettori et al. 2002). For our analysis we have 
taken a flat Universe with fl = 1.0, 17a = 0, and 
Hq = 50 kms~^ Mpc~^. In order to cover the whole 
cluster, we have used the observations of three satellites, 
namely Chandra, XMM-Newton and Beppo-SAX. Their 
views are complementary: Chandra's sharp angular res- 
olution (below 1 arcsec, corresponding to 1.7 kpc at the 
distance of A1795) is optimal for observing the core of the 
cluster. XMM-Newton has a broader PSF, but thanks to 
its wider field of view and large collecting area is partic- 
ularly well-suited for observing the intermediate regions. 
In the outskirts, where XMM-Newton^s view is contami- 
nated by a relatively high background and angular reso- 
lution is less critical, Beppo-SAX is a better choice. It is 
worth noting that the composite view from these satel- 
lites covers about two decades of the cluster extension, 
ranging from about 10 kpc to 1.5 Mpc. 

The 19.5 ks ACIS-S Chandra observation has been 
taken with a focal plane temperature of -110 degrees. 
The analysis has been performed with version V2.2 of 
the CIAO software and with the calibration data base 
CALDB 2.12. The details on the preparation and analy- 
sis of the XMM-Newton data are the same as in Molendi 
& Pizzolato (2001), the only difference being that here 



the new version V5.2 of the SAS software has been used. 
The preparation of the data from Beppo-SAX is de- 
scribed in De Grandi & Molendi (2001). 

In order to check the results of analytical deprojection, 
we have also undertaken the analysis with Ettori's matrix 
geometrical method. In the remainder of the paper we 
shall often address Ettori's method simply as "geomet- 
rical deprojection". The reader should remind, however, 
that this is only one of the possible geometrical depro- 
jection techniques. 

4.2. The geometrical deprojection 

We have selected 5 spectra from Chandra, 4 from 
XMM-Newton, 4 from Beppo-SAX, and have accumu- 
lated them in concentric rings centered about the X- 
ray emission peak. They have been analyzed using ver- 
sion 11.0.1 of XSPEC. A spectral analysis is common to 
the geometrical and the analytical deprojection. In the 
first case the analysis provides the 2-D emission- weighted 
quantities, still to be deconvolved. In the second, the 
spectral analysis directly yields the parameters of the 3- 
D distributions of the relevant physical quantities. 

Both spectral analyses include a wabs multiplica- 
tive component for the Galactic absorption (1.21 x 
10^° cm^^). A variable constant factor describes a mod- 
erate (< 10%) cross-normalization uncertainty among 
Chandra, XMM-Newton and Beppo-SAX. Since the spec- 
tra are not simultaneously fitted in the case of geomet- 
rical deprojection, the value of these cross-normalization 
factors must be provided. We have taken their values 
from the result of the (simultaneous) fit of the spectra 
performed with the smaug model described in the next 
Section. A further fixed geometrical constant factor has 
been introduced to account for the excision of the portion 
of the Beppo-SAX detectors shaded by the entrance win- 
dows. Since geometrical deprojection proceeds inwards, 
and necessarily only a finite number of rings can be used, 
the contamination from the emission of the cluster out- 
skirts must be considered. We have modeled t his c utoff 
effect as described by McLaughlin (1999). Table lGll sum- 
marizes the radii of the extraction regions of the spectra, 
their energy ranges, the geometrical corrections and the 
cross-normalization factors for the 13 datasets employed. 

The spectral analysis preliminary to the geometri- 
cal deprojection has been performed with a model 
wabs*constant*constant*mekal. The mekal component 
yields the emission integrals and the (projected) emission 
weighted temperature and metallicity of any ring, from 
which the relevant 3-D profiles can be extracted as de- 
scribed at the end of Section mi The error on these de- 
projected profiles is estimated with a Monte Carlo tech- 
nique. We build a sample of randomly distributed values 
around the mean, assuming Gaussian errors. From this 
sample a statistical collection of 1000 2-D profiles is ob- 
tained and deprojected: the standard deviation of the 
deprojected sample provides the required estimate of the 
error on the 3-D profiles. The outcomes relative to geo- 
metrical deproj ection a re s hown by the marks and error 
bars in Figures ICll IC2I and IC3I These Figures also plot 
the analytically deprojected results described in the fol- 
lowing Subsection. 

4.3. The analytical deprojection 
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In the language of XSPEC, our analytical model 
is wabs*constant*constant*smaug. The meaning of 
the first three components has already been discussed: 
smaug provides the parameters for the 3-D distributions 
of hydrogen density, temperature and metal abundance. 

The analytical model has 15 free parameters: two 
cross-normalization constants for the satellites and 13 
parameters for smaug. We have calculated the X-ray 
emission from A1795 with the Mekal emission code. The 
cutoff radius and the number of mesh points for the eval- 
uation of the integral have been fixed respectively 
to 2 Mpc and 15. We have checked a posteriori that 
the best-fit values of the parameters thus determined are 
not sensitive to these choices. In particular, changing 
the cutoff radius to 3 Mpc and rescaling the number of 
mesh points, the values of the best-fit parameters do not 
change appreciably, indicating that there are no impor- 
tant numerical effects associated with the evaluation of 
the integral Q. 

Not all the parameters available in smaug have been 
employed for the analysis of A1795: like all paramet- 
ric models, smaug calls for some a posteriori fine tun- 
ing of the parameters, suggested by the results of the 
fitting procedure. The exponent >f of the temperature 
gradient at intermediate radii during the fit systemati- 
cally pegs against its upper hard limit >c = 10, where we 
have chosen to freeze it. For this cluster the tempera- 
ture distribution does not depend much on the cooling 
radius rcooi, which is rather ill-determined by the fit pro- 
cedure. Therefore, we have decided to set it equal to 
the radius of the cluster density beta component rc. The 
polytropic index 7 of the plasma in the outskirts is poorly 
constrained by the Beppo-SAX datasets, whose statistics 
is rather poor with respect to Chandra and (above all) 
XMM-Newton. For this reason we have fixed its value, 
taking it from the statistical analysis by De Grandi & 
Molendi (2002) of the temperature profiles of a sample 
of Beppo-SAX clusters. These authors find an average 
7 — 1.2 for cool-core clusters. Accordingly, in our analy- 
sis we have fixed u)/ [5c = 7 — 1 = 0.2. 

The metal profile is a rapidly declining function of the 
radius, and therefore the parameter Zo is rather uncer- 
tain. When left free the radius rz assumes the unreliably 
small value of 1 kpc, comparable to the angular size of 
Chandra^s PSF at the distance of A1795. We have there- 
fore frozen rz to 8.24 kpc, i.e. one-half of the innermost 
Chandra^s radial bin. The parameter Zq assumes the 
physical meaning of the metallicity of the central bin. 

The deprojected profiles of density, temperature and 
metal abundance are plotted in Figures \CT\ IC2I and IC3I 
The best-fit values of the free parameters are shown in 
Table 1^21 with their 1-a errors. The reduced chi-square 
of the best fit is xVdof = 3652/3031 ~ 1.20, but its as- 
sociated probability is negligible, only 3.2 x 10"^^. This 
value is clearly underestimated: indeed, our analysis uses 
the spectral data from three different satellites, and thus 
a degree of cross-calibration error is unavoidable. We 
have partly kept notice of this by introducing two vari- 
able constants in our spectral analysis, but this is not 
good enough, since the spectral fit procedure tends to 
interpret the systematic cross-calibration differences as 
statistical errors, thus lowering the probability of the fit. 
We shall return in Section |31 to the problem of testing an 
analytical spectral model in presence of this effect. 



Like all models characterized by many fitting param- 
eters, also ours may be affected by the problem of sec- 
ondary minima. A clear example is the rather complex 
Ax^ pattern of the exponent Pc relative to the hydrogen 
density, plotted in the left panel of Figure lU^ The abso- 
lute minimum is well defined, but the plot is asymmetric, 
with several secondary minima. Complex behaviors are 
not always the rule: a much more well-behaved Ax^ is 
shown in the right panel of the same Figure, and refers 
to the temperature exponent ^. These examples demon- 
strate that it is always a good practice to check the re- 
sults against the possibility of secondary minima, before 
trusting blindly a deprojected profile. 



4.4. Inspection of the best fitting spectra 

We now turn our attention to the problem of checking 
the quality of our model for A1795. We have already re- 
marked that in our case the test is of relatively little 
use, being affected by the systematic cross-calibration 
effects arising from the use of data sets from different 
satellites. For this reason, in this and the next two Sub- 
sections we explore alternative routes. 

The first possibility, discussed in this Subsection, is es- 
sentially a qualitative examination of the residuals left 
over by the model on the spectrum of each annulus. Fig- 
ure IC5I plots the spectra of the 13 rings together with 
the related best fitting models and their residuals, ex- 
pressed as ratio data/model. As expected, the residuals 
of XMM-Newton^s data are smaller than the others: be- 
cause of its wide collecting area, XMM-Newton strongly 
dominates the statistics of the whole data set. 

The qualitative analysis presented here is easily under- 
stood recalling how the physical parameters shape the 
X-ray spectrum of a thermal plasma. The (squared) hy- 
drogen density provides the overall normalization of the 
spectrum, the plasma temperature fixes its exponential 
cut at high energies, and finally the metal abundance 
gives the equivalent width of a line. 

An inaccurate modclization of the profile would 
appear as a bad overall normalization of the spectrum 
of one or more annuli. The lower panels in Figure IC5I 
show such an effect in the innermost spectrum, indicat- 
ing that our model slightly underestimates the central 
density. The evidence of systematic normalization effects 
on the outer Beppo-SAX annuli is less clear. 

The systematic residuals at high energy are fairly small 
for any data set, showing that the 3-D temperature pro- 
file is modeled adequately by our functional form ^ . 

At the typical temperatures we have determined, the 
metal contribution mainly stems from iron. The local 
residuals about the 7 keV Fe-K line and the ~ 1 keV 
Fe-L blend test the accuracy of the determination of the 
metal distribution. The Fe-K line has a good statistics 
only where the cluster is hot enough, i.e. in correspon- 
dence of the area covered by XMM-Newton. Near the 
core, where the cluster is cooler, the line emission is dom- 
inated by the Fe-L blend. Unfortunately, the Fe-L blend 
statistics is not high, and Chandra^s determinations of 
the me tal abundance are rather scattered, as shown in 
Figure IC3I In the absence of strong constraints from 
the data, therefore, we conclude that even for the metal 
distribution our model is fairly reliable. 
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4.5. Comparing the analytical and geometrical routes: 
a discussion 

A second way to check the goodness of our fit is to com- 
pare its resuhs with those from the matrix geometrical 
deprojection. 

The best-fit analytically deprojected 3-D profiles of hy- 
drogen density, te mperature a nd metal abundance are 
plotted in Figures ICll IC2I IC3I together with their geo- 
metrically derived counterparts. Each of these Figures 
also shows the 1-cr error band about the best-fit profile. 
The band enclosing the best-fit curve is the area spanned 
by all the functions whose parameters yield a difference 
< 1 with respect to the best-fit x^- This method 
of characterizing the errors is thoroughly described in 
Appendi x C. 

Figure [HTI shows that our modeled density is smaller 
than the geometrically deprojected value in the central 
annulus, and larger in the two outermost rings. Else- 
where the agreement between the two techniques is good. 
These results agree with the qualitative analysis of the 
model spectral residuals outlined in the last Subsection. 
Leaving aside the outermost Beppo-SAX annuli, whose 
statistics may be weak, the combined results of analyti- 
cal and geometrical deprojection seem to suggest that - 
at least for A1795- the density profile is somewhat more 
peaked than a double beta-model. 

ft is apparent from Figure[n2]that the 3-D temperature 
suffers a sharp discontinuity at about 37 kpc. This fact 
was anticipated by the strange behavior of the parameter 
X noticed in Section ft is also evident from the geo- 
metrical deprojection, and was already noticed by Ettori 
et al. (2002) in the same position. We view the abil- 
ity of our model to spot such an unexpected distribution 
in the three-dimensional temperature profile positively. 
Our deprojection finds a maximum 3-D temperature of 
5.8 keV, to be compared with 7.2 keV of the geometrical 
deprojection. This difference is due to discr epancies in 
the calibration of the satellites. Figure IC6I shows that 
the temperatures measured by XMM-Newton seem too 
low to match well the neighboring values provided by 
Chandra and Beppo-SAX. Because of their strong statis- 
tical weight, XMM-Newton^s measures tend to lower our 
deprojected profile, reducing also the maximum temper- 
ature. 

Figure E^l compares the 3-D metal profiles. The agree- 
ment between the analytical and the geometrical depro- 
jection is good. In this case it is particularly clear how 
our profile turns out to be smoother than its wiggly geo- 
metrically deprojected counterpart, which would need a 
regularization in view of deriving other physical quanti- 
ties. 



4.6. Back to Flatland 

The third reliability test we suggest is very stringent, 
since it compares directly the model predictions with the 
emission weighted, observed quantities. It consists in re- 
projecting the 3-D best-fit model profiles and comparing 
them to the 2-D measurements. 

The procedure for temperature T and metal abundance 
Z are the same, so we sketch it for temperature only. 
The temperature T2D{bi,b2) averaged over the volume 



seen under the annulus bounded by the rims bi and 62 is 



T2D{blM) 



j'^ldbj^dlj{r)nr) 

j,';dbj,-dij{r) - 



where j(r) is the plasma emissivity. The integrals in dl 
along the line of sight can be expressed as radial integrals 
with the geometrical relation r'^ — P + The resulting 
double integrals (in db and dr) can be reduced to sin- 
gle integrals, as described in Appendix A. The previous 
formula is therefore 



XT dr /C(r, bi, 62) j{r) T{r) 
drlC{r,bi,b2)j{r) 



(10) 



where the kernel K. is given by Equation ©. As long as 
the temperature of the ICM remains above ~ 2 keV, the 
contribution from line emission is negligible and the emis- 
sivity is mainly due to thermal Bremsstrahlung, yielding 
j (X n^^T^/^. We calculate the emission weighted quan- 
tities under this assumption. 

Figure [UgI plots the calculated 7215(61,62) binned over 
rings i7 kpc wide over the observed values, directly ob- 
tained from the mekal analysis. Similarly, Figure IC7I 
compares the measured and re-projected metal abun- 
dance profiles. In both cases, the results are in good 
agreement with the observations. 

For the re-projection of hydrogen density we have 
taken a slightly different route. The geometrical depro- 
jection spectral analysis yields the mekal normalization 
of the spectrum of each ring (Arnaud & Dorman 2000) 



M 



10^ 



' J dVUf-riH 



where the integral is extended over the volume seen 
through the ring. Da and z are respectively the angular 
distance and the redshift of the cluster. From our analyt- 
ical deprojection the distributions nn and Ue — 1.21nH 
of hydrogen and electron densities are known, and there- 
fore the quantity M can be evaluated for each ring, and 
compared to the corresponding value provided by the 
mekal analysis. Figure [HHI shows the results, which again 
are in good agreement, with the possible exception of the 
innermost ring. 

5. ON THE REJECTION OF AN UNSATISFACTORY MODEL 

So far we have devoted a good deal of attention to 
check the goodness of a single fit. In this Section we 
explore the related issue of how to check the goodness of 
a fitting function. The problem is to assess our ability to 
reject an unsatisfactory spatial profile. We would like to 
consider only self-contained tools, which do not require 
to compare our results with those from other kinds of 
deprojection. 

It is clear that the spectral test and the reprojection 
test discussed above are well-suited for this goal. A spa- 
tial functional form ought to be rejected when it leaves 
relevant systematic errors on many annular spectra and 
on the re-projected quantities, and when these residuals 
cannot be reduced with any choice of the functions pa- 
rameters. For instance, the reprojection test rules out 
an isothermal profile for Ai795: the projection of a con- 
stant function is again a constant function, clearly at 
odds with the 2-D profile shown in Figure lUHl The same 



7 



result comes from the spectral test, since the systematic 
residuals left over by an isothermal profile are large. 

Both the reprojection and the spectral test have been 
devised to complement the test when systematic 
cross-calibration effects due to the use of inhomogeneous 
data sets are important. Obviously, these effects lessen 
the reliability of the test also in testing the quality 
of a fitting function. There is an important limit, how- 
ever, in which the test is still extremely useful, namely 
when the function is grossly wrong. In this case, the sta- 
tistical errors yield a much larger contribution than the 
systematics, and the test is reliable. For instance, had 
we taken a simple power-law for the nn distribution, we 
would have obtained xVdof = 62271/3035 ~ 20.5. The 
badness of the fit is so marked that the value of the x^ 
is only marginally affected by systematic errors, and can 
be confidently used to rule out such an unlikely profile. 

The x^ test also rules out both an isothermal, a beta- 
model and a polytropic profile for the temperature 
The reduced chi-squared is xVdof = 4052/3034 ~ 1.3 
for the isothermal model, xVdof = 3740.6/3031 ~ 1.2 
for the beta-model and xVdof = 3781.1/3034 - 1.24 for 
the polytropic model. These figures are not exceedingly 
larger than our best-fit value, nevertheless the models 
can be s afely rejected with an F-test, summarized in Ta- 
ble (|C3|) . As a by-product, these tests indicate that mod- 
els with few parameters cannot reproduce adequately the 
temperature profile of A1795, and all the free parameters 
employed to describe our best-fit temperature profile are 
necessary. 

We have also tested the modified double beta-model 

riH = no {/[I + {r/m)T^- + (1 - /)[1 - {r/rcff^}. 

(11) 

put forward by Ettori (2000) to correct the original dou- 
ble beta-model, which lacks a satisfactory physical mean- 
ing. 

The rationale of this model is the confinement of the 
cool gas within the sphere of radius rc where radiative 
cooling is important. We have tried several combinations 
of parameters: letting all them free, linking or freezing 
some of them, but in any case the corrected double beta- 
model has somewhat worsened the quality of the fit with 
respect to the native double beta-model. For instance, by 
freeing all the density, temperature and metal abundance 
parameters (with the exception of the metallicity radius 
rz and the temperature exponent >c, which is irrelevant 
to the present discussion), the x^ has worsened from 
xVdof = 3652.3/3031 ~ 1.20 (original double-beta) , to 
xVdof = 3785.1/3030 ~ 1.25 (corrected double-beta). 
The best-fit confinement radius rc — 180 kpc is essen- 
tially inside the XMM-Newton and Beppo-SAX datasets, 
which are therefore described by the hot-beta compo- 
nent alone (the first right-hand term in Equation lllp. 
In other words, we have given up the degrees of freedom 
of the cold component to model the most statistically 
relevant datasets, which of course reflects on the final 

The use of the x^ test to reject a model gets difficult 

The beta-model for the temperature has been obtained by fix- 
ing some parameters in the general profile JSJ , and corresponds to 
the function T = Tofc(r) in the notation of that Equation. The 
polytropic profile T oc has been obtained with a slight modifi- 
cation of the smaug code. 



as the model improves, i.e. the functional form is flexible 
enough to be able to reproduce fairly well the qualitative 
behavior of the actual 3-D distributions. In this case the 
role of statistical errors decreases in front of the system- 
atics, which tend to "saturate" the x^j especially when 
dealing with inhomogeneous data sets. In order to get 
rid of this effect and focus on testing the reliability of our 
functional forms, we have analyzed separately Chandra 
and XMM-Newton datasets. In both cases the reduced 
X^ diminishes, (xVdof = 1362.5/1098 ~ 1.24 for the 
5 Chandra datasets and xVdof = 1806.2/1535 ~ 1.18 
for the 4 XMM-Newton rings), but the associated prob- 
abilities remains small. Since calibration effects are re- 
duced, they are unlikely to be fully responsible for this. 
This indicates that the introduction of a functional form 
is certainly an oversimplification of the complexity of a 
system like a galaxy cluster. It smoothes out local in- 
homogeneities, deviations from spherical symmetry, and 
so on. It is clear that if the spectra have good statistics 
as in the present case, these unmodeled effects refiect on 
the final value of x^- 

We conclude that a bad x^ model is not always to be 
discarded. If it has good reprojection and spectral resid- 
uals tests, it should be regarded as a reasonable model 
for the 3-D structure of a cluster, and may be as useful as 
a profile obtained via geometrical deprojection. Indeed, 
the deprojected profiles are usually taken as a starting 
point for other calculations, like the gravitating mass of 
the cluster. In a geometrical deprojection one may ob- 
tain a good x^ statistics on the spectral fits, but the 3-D 
raw proflle often needs smoothing to be used for further 
analyses. It is clear that smoothing erases a posteriori 
some of the original information. Along the analytical 
route, some of this information is disregarded a priori, 
at the cost of a relatively high x^ • The 3-D proflles, how- 
ever, are smooth and ready for further uses. Moreover, 
the knowledge of the functional forms of the distribution 
in some occasion may prove useful. 

To summarize, some warnings are in order: the analyt- 
ical deprojection is not foolproof: a preliminary analysis 
of the regions to deproject is necessary to be sure not 
to miss important details of the structure of the cluster 
to be analyzed. Moreover, the flnal fit should never be 
taken as gospel truth, but always checked against pos- 
sible flaws in the model or in the functional forms for 
the spatial distributions of the physical quantities. To 
this end, the x^ test of the results should be used with 
care, considering the possible role of systematic cross- 
calibration effects. Generally speaking, our advice is to 
privilege the spectral and the reprojection tests to check 
the quality of a flt, as well as of the functional form of a 
fitting function. 

In conclusion, if the analytical method is used carefully, 
it may provide a reliable deprojection of a galaxy clus- 
ter. Further, the knowledge of explicit functional forms 
for the distributions of density, temperature and metals 
somewhat simplifies the after-deprojection task of work- 
ing out indirect quantities like the mass profiles of the 
cluster, as explained in the next Section. 

6. RECOVERING THE MASS PROFILES 

In this Section we show how the results obtained from 
analytical deprojection can be used to deduce the mass 
profiles of the cluster. This part of the procedure is not 
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implemented directly in the model smaug within XSPEC, 
but it represents an important complement to the infor- 
mation provided by the spectral deprojection. Indeed, 
it is a sort of a posteriori justification of the whole ma- 
chinery of analytical deprojection. Therefore, the codes 
necessary to this part of the spatial analysis of a clus- 
ter will be made available to the community to complete 
smaug. 

Once the profiles of hydrogen density, temperature and 
metal abundance have been determined, in principle it is 
easy to recover other interesting physical quantities like 
the gravitating mass, the X-ray emitting gas mass frac- 
tion or the metal mass. This is an advantage of our ana- 
lytical deprojection: the knowledge a priori of the func- 
tional forms of uh , T or Z allows an easy automatization 
as well as the numerical stability of all the necessary pro- 
cedures. The hypothesis of hydrostatic equilibrium and 
the equation of state for the ICM yield the profile of the 
gravitating mass within radius r: 



Mg{< r) 



kTir) fdlogUH d logT 



Gfinip \ d log r d log r 



(12) 



(see e.g. Sarazin 1988). The quantities on the right-hand 
side of Equation H12() can be calculated analytically from 
the profiles © and © : some lengthy but straightforward 
algebra yields 



d log UH 

d log r 



d logT 



/tIc + (1 - /)% 



1 + {rjry 



■ ric + 



(1 - m 

l + (r,/r)2"% 



-2uj 



dlogr l + {rt/r)^ To + Til{r) 



arctan \(jlr\^of \ 1 + (r/riso)^'* 1 + (r/rcooi)^ 

In the case of A1795, as we have anticipated, the hy- 
pothesis of hydrostatic equilibrium breaks down at about 
37 kpc from the center, because of the discontinuity 
in the temperature profile. Mathematically, in Equa- 
tion (|12|) the derivative of T is large and positive, thus 
making Mg{< r) negative. This is shown in Figure lUOl 
by the plots of the derivative of T and the calculated 
Mg{< r) profile. It is important to remark that the hy- 
drostatic determination of the gravitating mass only fails 
near the temperature discontinuity. The physics of this 
region is complex: according to Markevitch, Vikhlinin & 
Mazzotta (2001) the cold gas is slowly moving within the 
gravitational potential of the bulk of the cluster, proba- 
bly after a past subcluster infall. Elsewhere the hydro- 
static condition holds, and this recipe for the determina- 
tion the gravitational mass is relatively safe. We remark 
that the ability of our model to spot an unexpected be- 
havior of the gravitating mass ought to be considered an 
index of its flexibility, rather than a drawback. A simple 
numerical integration of the hydrogen density pro- 
vides the mass profile of the X-ray emitting gas enclosed 
within the radius r: 

Mx{<r) -1.39 nip [ 4tt dr' r'^nH{r'), (13) 
Jo 

where the factor 1.39 is the correction from the hydro- 
gen to the average ion density. This mass profile is de- 
fined by an integral, thus it is always well defined, unlike 



the gravitating mass, whose determination relies on the 
derivatives of potentially bad- behaved functions, as hap- 
pens for A1795. Figure ICIOI shows the X-ray emitting 
gas fraction /(< r) — Mx{< r)/Mg{< r). It reaches 
a plateau of about 0.17, but slightly rises in the out- 
skirts. This final trend (which is nevertheless consistent 
at the 2-cr level with being a constant) is related to the 
possible overestimate of the hydrogen density profile at 
large radii previously noticed. This leads to overestimate 
the gas mass profile Mx{< r), and hence the fraction 
fi<r)=Mx{< r)/Mg{<r) . 

The combination of the hydrogen and metallicity pro- 
files provides the metal distribution. At the temperatures 
observed in A1795, the main contribution to the metal- 
licity profile Z{r) comes from iron, the concentration of 
which is approximately 

nFo(r) ~ [Fe/H]nH(r)Z(r), (14) 

where [Fe/H] — 4.68 x 10^^ is the solar abundance of 
iron with respect to hydrogen, according to Anders & 
Grevesse (see Arnaud & Dorman 2000 and references 
therein). From riFe it is easy to deduce the enclosed iron 
mass Mfc(< r). Figure lUTTI plots these quantities. 

7. SUMMARY 

In this paper we have described a new technique for 
recovering the three-dimensional structure of clusters of 
galaxies. The spectrum of a cluster is modeled as a func- 
on of a number of parameters, whose values are de- 
termined by fitting the model to the observed spectrum. 
This model can be made to work sufficiently fast as a sub- 
routine of XSPEC because the calculation of the spectral 
model calls for the evaluation of a single rather than a 
double integral (Sect ion [3 and Appendix A), and because 
a further analysis has shown that a time-consuming eval- 
uation of the spectral emission function at several points 
throughout the cluster is unnecessary (Appendix B). 

The parameters of the model define the known ana- 
lytical functions of the spatial distributions of hydrogen 
density, metal abundance and temperature. From these 
quantities, under suitable assumptions it is straightfor- 
ward to obtain other physical quantities like the gravi- 
tating mass, the X-ray emitting gas fraction or the dis- 
tribution of metals within the cluster. 

The distinctive feature of this technique is the possi- 
bility of adopting known analytical fitting functions for 
the use of deprojected data. After the spectral analysis, 
no further work is required to obtain and regularize the 
three-dimensional profiles. Noticeably, our method has a 
relatively simple way to characterize the uncertainty of 
the results (Appendix C), and does not require a Monte 
Carlo technique. The use of known functions simplifies 
the reconstruction of the cluster 3-D structure after the 
spectral deprojection. Once the best-fit parameters have 
been worked out from the spectra -with all the caveats 
outlined in Section [S]- this part of the procedure can be 
easily and safely automatized. 

We have shown that the functional forms of nn , T and 
Z, if suitably chosen, may encompass a wide range of pos- 
sible behaviors of the profiles of these physical quantities. 
This was shown by the ability of our model to spot a dis- 
continuity in the temperature profile of A1795. If the 
data from the last generation of X-ray satellites required 
it, however, it would be extremely simple to modify these 
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functional forms to achieve a closer agreement with new nique may provide a valid and useful alternative tool for 

observations. the deprojection of galaxy clusters. 

Finally, we have shown that the test, the reprojec- 
tion test and the analysis of the spectral residuals left 

by the models are useful tools to check the reliability The authors wish to thank Keith Arnaud and Ben 

of the results, and possibly may suggest new and bet- Dorman at HEASARC for their prompt and kind guide 

ter functional forms for the 3-D distributions of density, through the mysteries of XSPEC. Thanks also to our ref- 

temperature and metals. eree Stefano Ettori, whose suggestions greatly improved 

In conclusion, we believe that our deprojection tech- the original manuscript. 

APPENDIX 

REDUCTION OF THE DOUBLE INTEGRAL 

In this Appendix we derive Equation Q of the text, starting from Equation The integral therein is extended 
over the solid angle subtended by an annulus bounded by the rims 6min and 5niax- First make the change of variable 
from the solid angle to the projected radius <Kl = 2 7Tbdb/D\, where Da is the angular distance of the source. Then 
substitute Equation (O into QJ, taking into account the invariance law I'^,/E'^ = Ie/E^ for the specific intensity. 
The result is the photon flux seen by the observer 

47r 

where the quantity ce' = jE'/E' is the photon emissivity at the source, in photons cm^'^ s^^ keV^^. The primes 
indicate that th e en ergy with respect to the source and the observer's rest frame differ for the redshift factor \ + z. 
The integral in I|A1|1 can actually be reduced to a single integral. By inserting a couple of theta functions 

'0 for a; < 
,1 for a; > 

the double integral reads 



= TW7T-7^ I dr^^^ ^E'ir), (Al) 



(x) = 



dbb dr eE-{T)^==e{r - b) e{r - 6„,in), 

f>min "'0 



r 

where the order of integration can be reversed; the integral in db can then be explicitly solved, yielding 

1 ( (r"^ - b ■ ^"1^/^ for r < 6 

77 7 ^ u\ al J, ^ J V "mm j ^'^^ ' ^ "max 

On substituting this in the integral over dr, formula Q of the text is recovered, QED. 

NUMERICAL EVALUATION OF THE SPECTRUM 

This Appendix provides some technical details about the numerical calculation of the integral in Equation 10} of the 
text. It is useful to single out the hydrogen density in the emissivity function: 

JE Ae{T, Z) _ 2. (rj. 

€E = — = rieUH — = njjXE{T, Z). 

Omitting irrelevant factors, the integral to be calculated is 

= dr^XEiT,Z), (Bl) 

where the differential emission integral 

^ = JC{r; 6max) njj{r), 
dr 

is far more sensitive to the radius r than the spectral function \e{T, Z), which depends implicitly on r through the 
functions T(r) and Z{r). This suggests to split the integration interval in several isothermal shells with homogeneous 
metallicity. Therefore, in each shell the function Xe{T,Z) needs only one evaluation. A high numerical accuracy is 
put only in evaluating the emission integral of the shell. This approach greatly reduces the (computationally heavy) 
evaluation of \e{T, Z) in several mesh points, but preserves the accuracy of the results. 

The range of integration is cut-off with the introduction of a user-defined maximum radius Tmax, typically in the 
order of ~ 2 Mpc, and the interval [6min, ?'max] is spht into N shells bounded by the nodes rp , ri , ■ ■ • .rjv (in general 
not evenly spaced); the extremes tq and r^r coincide with 6min and Tmax respectively. Equation (IB1|1 is rewritten as 

(B2) 



r dr 



dr- 
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where the average spectral function Xsii) over the i-th emitting shell is defined as 



dr 



(0 



(B3) 



Up to now the equation are exact: the approximations are to be introduced now in the evaluation of the average 
spectral function A^. For the sake of readability of the formulae we drop both the spectral index E and the label of 
the volumes. In Equation IIB3I) we expand A(r, Z) in a neighborhood of the (still undefined) average temperature and 
metallicity T and Z, defined on each shell: 



A(r, z) = A(r, z) + |A(r - T) + ^{z z) 



The term denotes higher-order terms in (T — T) and {Z — Z). If we put this expansion in the right-hand side of 
llB3ll we obtain 



A 



fdr^T 

J dr 



T 



aA 

dZ 



Z 



If we define the shell averages T and Z as 



T : 



_ f dr^ Z 

^ — J dr 



(B4) 

(B5) 
(B6) 



we find A = A(T, Z) + , i.e. the (approximated) spectral function evaluated in T and Z is first-order accurate with 

respect to the (exact) spectral function (|B3p . emission- weighted within a shell. 

The trapezoidal rule is accurate enough for the calculation of the integrals in the definitions ljB5(l and l)B6|l : on the 
shell [ri,ri_i], for instance 



T. 



IC{i)n%{i) + IC{i-l)njj{i-l) 



where IC(i), T(i), and nnii) denote the values of the functions JC, T, and nn evaluated at r^. As mentioned before, 
the evaluation of the emission integral 0(i) in the sum (jB2p calls for more accuracy, and therefore we adopt a 10-point 
Gauss-Legendre formula (see e.g. Press et al. 1992). In each integration shell the differential emission integral is 
calculated 10 times, which is accurate and relatively fast; the evaluation of the spectral function Xe{T,Z), which is 
substantially slower, occurs only once. 

A final remark is in order about the spacing of the integration nodes ro,ri, • • • .r^. Since the differential emission 
integral usually strongly peaks at r < 6max, we have chosen a quadratically spaced grid, which yields a better accuracy 
than an evenly spaced one. 



ON THE CHARACTERIZATION OF THE ERRORS 

In this Appendix we sketch our method to characterize the error on a function /(r; a) dependent on the radial 
coordinate r as well as on the iV parameters a = (oi, • • • ,ai\r). Our problem is to show how their uncertainties reflect 
on the function /. We start by considering the Taylor expansion of the about the best-fit values of the parameters 
a ~ (fli, • • • , ojv). Here has a minimum, and its expansion to the lowest order reads 

A,».,»(a)-X^(a) = iE(''-^).('-^)'^- («) 

Z,fc = l * 

In the same fashion, the first-order Taylor expansion of /(r; a) about a reads 

A/^/(r;a)-/(r;a) = f](a-5),^:^^. (C2) 

i—l ^ 

In what follows we assume that r is fixed, and focus our attention on the iV-dimensional parameter space. Her e the 
manifold Ax'^ = constant is the {N — 1) - dimensional surface Ti. of an ellipsoid. In the same space, Equation ljC2|l 
defines a family V of parallel planes dependent on the intercept A/. The intersection set HnV contains the parameters 
which for the assigned value of Ax^ give the absolute difference | A/| with respect to the best-fit determination of /. If 
I A/ 1 is large, the intersection is empty: the error on / is too large to be consistent with the given Ax^. On the other 
hand, if |A/| is small, there are infinitely many points in the intersection, i.e. several combinations of parameters are 
able to define many functions characterized by the same Ax^. It should be clear (and may be proved rigorously) that 
the largest value of |A/| compatible with the fixed value of Ax^ occurs for the values q of the parameters for which 
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the ellipsoid (jCip and one plane of the family (|C2|I a re tangent. Quite generally, given a point T of coordinates q lying 
on the surface Ax^ = constant of the ellipsoid ljCl|l . the plane tangent to the ellipsoid at T is 

N 

Y,{E-%k{q~~aUa-~a)u = ^x\ (C3) 

i,k=l 

where the error matrix £i k is the inverse of the matrix defined by 

' 2 daidak 

The tangent plane ljC3p at T coincides with the plane of the family ljC2|l passing through T if the coefficients of both 
are proportional to each other. This condition yields 

A/= ^ t ..^ILJL^ (C4, 



Since £ is the reciprocal of a (symmetrical) Hessian matrix evaluated at a local (and hopefully, absolute) minimum 
of x^, it is symmetrical and positive definite; the equation for A/ has two roots ±|A/|. As the radial coordinate r 
varies, /(r) ± A/(r) dra ws th e error band associated with the given value of A^^- 

Alternatively, formula ljC4p can be obtained by maximizing the value of |A/| (as a function of the parameters a), 
subjected to the condition that a belongs to the assigne d m anifold Ax^ — constant. An application of the standard 
Lagrange multipliers technique leads again to Equation ljC4|) . 

In Equation (|C4|I the components of Va/(?'; a) are known from the functional form of / and from the best-fit values 
of the parameters. The matrix £ is not directly known, but it can be worked out using the variances and the principal 
axes returned by XSPEC at the end of the fit procedure. The orthogonal matrix V of the principal axes is made up by 
the (row) eigenvectors of f , and the variances are the corresponding eigenvalues. Therefore, £ is given hy £ = *VI]V, 
where S is the matrix with the variances on its main diagonal, and *V is the transpose matrix of V (the interested 
reader is referred to any textbook of Hnear algebra, e.g. Lang 1966). 

The next degree of complication occurs when more than one function is present, as it is our case (the reader is 
reminded that we determine simultaneously the functions n//, T and Z). Suppose that besides / there is another 
function g: for instance / and g could be the temperature and the hydrogen density. In order to avoid unnecessary 
complications, suppose that each function depends only on the radius r and on one parameter: / = /(r, a) and 
g = g{r,b). The parameter space is the two-dimensional plane ab, the error matrix has rank 2, the x^ depends both 
on a and on b and the surface Ax^ = constant is an ellipse in the plane ab. Assume for the moment that a and b are 
mutually independent, i.e. they are not linked to each other in the fitting procedure: notice however that this by no 
means entails they are statistically uncorrelated. As the parameter a varies about its best fit value a, the error A/ 
also varies, and in the same fashion as before its largest value 

Af = Ax'£^^(^X (C5) 



occurs where the error ellipse and the line 

A/ 

are tangent. Similarly, the error on g is 



5/ 
da 



da 
(a — a) 



Ag^^Ax'£2 2\^^j . (C6) 

It might be objected that this procedure does not take into proper account a possible statistical correlation between a 
and 6. This is expressed by the off-diagonal terms of the error matrix, that however do not appear in formulae lj(^5|l and 
(|C6ll . This objection is incorrect, since the correlation between a and b (if any) intervenes also through the diagonal 
terms £ii and £22- Indeed, the correlation introduces a tilt in the principal axes of the error ellipse with respect to 
the c oordinate axes of the parameter plane. Therefore the diagonal terms of the error matrix in Equations (|C5|) and 
(fC6l differ in the cases of correlation or uncorrelation between the fit parameters a and b. We may conclude that the 
errors Af and Ag are correctly calculated in both cases. 

The situation is slightly more complicated if an explicit link b = 4>{a) is introduced in the fit between a and 6; b is 
no more a free fit parameter, and its role is taken over by a. Because of the link, the quantity dg/da does not vanish, 
since 

^9 9g , 
d-a=db'^^''^- 

The formula (|C6p in this case is replaced by 

Ag"" ^Ax^ \£ii[^<P'] +2£i2^0'+f22 ( ^ 1 !>. (C8) 
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This procedure of dealing with the hnks between parameters aUows to derive the error for one special parameter, 
namely the hydrogen central density no of the model smaug. no is special because during the fit it is formally fixed, 
but its role is played by the normalization J\f of smaug, being Af = Uq. We interpret this relation as a formal link 
uq — 4>{M) — \/M between TV and no. In Equation (|C4|) . the derivative with respect to the fit parameter A/" of a 
function / dependent on no reads 

As a final point, let us consider how it is possible to characterize the error of a function h = h[f{r, a), g{r, h)) dependent 
both on / and g: for instance, if / and g are the temperature and the hydrogen density h — f g is the pressure. The 
basic formula is always (|C4|I . and therefore it is only necessary to evaluate correctly the derivatives on its right hand 
side. If the fit parameters a and b are not linked, it is dh/da = {dh/df ){df /da) and dh/db = {dh/dg){dg/db). The 
occurrence of a link b = (j>{a) leads to 

dh dhdf dhdg . 
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lO'O g_OT ^_0T FT T 6'0 

[ rao] oi^iej 

Fig. CI. — Upper panel: the deprojected hydrogen density profile. The central continuous curve is the best-fit profile yielded by our 
analytical deprojection with the double beta-function (see Equation The other curves are the upper and lower boundaries corresponding 
to a one-sigma error in the parameters (see text). The markers show the data and error bars from Ettori's matrix geometrical deprojection. 
Circles correspond to Chandra data, triangles to XMM-Newton and diamonds to Beppo-SAX. Lower panel : the ratio between the hydrogen 
density derived from the geometrical and the analytical deprojection. 
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Fig. C2. — Upper panel: the depr ojected temperature profile, whose functional form is given by Equation JHJ. The meaning of the 
symbols is the same as in Figure lUll Lower panel : the ratio between the temperature derived from the geometrical and the analytical 
deprojection. 
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[s^Tun ji9|os] aou^punqig I^'ISH OT^ej 

Fig. C3. — U pp er p anel: the deprojected metal abundance profile, with functional form j5J. The meaning of the symbols is the same as 
in Figures ICll and )C2l Lower panel : the ratio between the metal abundance derived from the geometrical and the analytical deprojection. 
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Fig. C4. — The Ax^ on two fit parameters, obtained by running the XSP-BC coniniand steppar. Left panel: the Ax^ on the exponent /3c 
of the hydrogen cluster component. Despite the presence of secondary minima, the best-fit value of this parameter is robust. Right panel: 
the same plot for the exponent { of the temperature profile shows a much smoother behavior. 



Table CI. Extraction regions for the spectra 



Instrument 


Cross 
normalization 


Spectral range 
[keV] 


Inner ring 
[arcsec] 


Outer ring 
[arcsec] 


Shading 

factor 


Chandra 


106% 


0.5-7.0 
0.5-7.0 
1.0-7.0 
1.0-7.0 
1.0-7.0 



10 
25 
40 
60 


10 
25 
40 
60 
90 


1.000 
1.000 
1.000 
1.000 
1.000 


XMM-Newton 


100% 


1.0-10.0 
1.0-10.0 
1.0-10.0 
1.0-10.0 


60 
120 
180 
240 


120 
180 
240 
300 


1.000 
1.000 
1.000 
0.981 


Beppo-SAX 


95% 


2.0-10.0 
2.0-10.0 

3.5 10,0 
2.0 8.0 


240 

360 

480 
720 


360 
480 

720 
900 


1.000 

0.981 

1.000 
0.791 



Note. — The cross-normalization constants among the instruments are free parameters of the model. The value relative to 
XMM-Newton is taken as unity by definition. The shading factor takes into account the reduced area of the extraction region 
of two Beppo-SAX spectra due to the entrance window. 



Table C2. Results for the smaug model. 



Free parameter 


Best-fit value 


One-sigma error 


Units 


/ 


4.41 X 10"^ 


1.40 X 10-2 




no 


2.89 X 10-2 


8.91 X 10-4 


cm-3 


/3c 


1.76 


7.26 X 10-1 




rc 


5.03 X 10-1 


1.58 X 10-1 


Mpc 


Pa 


8.76 X 10-1 


5.72 X 10-2 




rg 


5.07 X 10-2 


2.40 X 10-^ 


Mpc 


Zo 


1.35 


1.34 X 10-1 


solar 


C 


2.05 X 10-1 


1.60 X 10-2 




To 


2.85 


1.27 X IQ-i 


keV 


Ti 


6.55 


6.20 X 10-1 


keV 




3.75 X 10-2 


3.37 X 10-3 


Mpc 




3.94 X 10-1 


9.13 X 10^2 




nail 


1.26 


5.80 X 10-1 


Mpc 



Note. — The parameters not listed are the column density of Galactic absorption uh = 1.21 x 10 cm ^, x frozen to 10, 
rz frozen to 8.24 kpc i.e. one-half of the innermost radial bin, r-cooi (set equal to rc) and w = 0.2 /?c. 
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0" < b < 10" 



10" < b < 25" 



25" < b < 40" 




0.5 1 



40' < b < 60" 



channel energy (keV) 




0,5 1 



channel energy (keV) 



60" < b < 90" 



channel energy (keV) 




160" < b < 240" 



channel energy (keV) 



360" < b < 4B0" 



channel energy (keV) 



channel energy (keV) 



' < b < 130' 



channel energy (keV) 




460" < b < 720" 



channel energy (keV) 



channel energy (keV) 



Fig. C5. — For any ring, the best-fit spectrum plotted over the data (upper panels). The lower panels show the residuals of the model, 
plotted as ratio data/model. Panels 1-5 refer to Chandra, 6-9 to XMM-Newton and 10-13 to Beppo-SAX. 
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Fig. C6. — Upper panel: the re-projected temperature profile plotted over the corresponding 2-D measures provided by the mekal 
analysis. Lower panel: the ratio between the measured emission-weighted temperature and the corresponding rc-projccted values. 



Table C3. Checking the temperature distribution. 



Profile 


Free parameters 




Degrees of freedom 


F-test statistic value 


isothermal 


To 


4050.8 


3035 


82.7 


beta law 


To nail a; 


3731.1 


3033 


32.7 


polytropic 


t 


3781.1 


3034 


35.6 


best-fit 


To Ti nso £, rtail 


3652.3 


3031 





Note. — The number of free parameters required to describe correctly the temperature distribution. In the beta-model the 
exponent is free. The F-test statistic values compare the model of each line with respect to the best-fit one. In the row of 
the polytropic temperature model the names of the free parameters have not been reported, since the function ||SJ| has been 
modified. In this case, the temperature 3-D distribution is described by two parameters. In all cases the associated probability 
is well below 10"^'', and has not been reported. 
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[s-^iun jT3];os] Aipni^T-^lAt otitsj 

Fig. C7. — Upper panel: the re-projected metallicity profile plotted over the corresponding 2-D measures. Lower panel: the ratio between 
the measured emission-weighted metal abundance and the corresponding re-projected values. 
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Fig. C8. — Upper panel: The re-projected profile of the emission integral EI = J dVuenn plotted over its corresponding 2-D measured 
counterpart. In the y axis Dy^ = 3.4 X 10^ Mpc and z = 0.0631 are respectively the angular distance and the redshift of A1795. As before, 
circles indicates Chandra data, triangles XMM-Newton data and diamonds Beppo-SAX data. The observed quantities stem directly from 
the normalization of the spectra yielded by the mekal analysis. The re-projected quantities have been calculated by evaluating numerically 
the emission integral from the best-fit hydrogen density profile. Lower panel: The ratio between the mekal emission integral and our 
reprojected value. 
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Fig. C9. — The gravitating mass calculated with Equation 1121 (upper panel). Because of the discontinuity edge in the temperature 
profile, the hypothesis of hydrostatic equilibrium breaks down at ~ 37 kpc. This is best appreciated in the lower panel where we report 
the logarithmic derivative of the temperature profile. 
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Fig. CIO. — The X-ray emitting gas fraction within radius r is the ratio of the masses given by 1131 and 1121 . The plot excludes the 
innermost region, where the hydrostatic equilibrium fails. The slight rise in the outskirts is statistically significant only at the l-a level, as 
discussed in Section ISl 
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Fig. Cll. — The iron density (upper panel) and the enclosed iron mass within radius r (lower panel). The error curves correspond to 
l-O" uncertainty in the parameters, and have been calculated as explained in the text (Appendix C). 



